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We report experiments on the full space and time resolved statistics of capillary wave turbulence 
at the air- water interface. The three-dimensional shape of the free interface is measured as a function 
of time by using the optical method of Diffusing Light Photography associated with a fast camera. 
Linear and nonlinear dispersion relations are extracted from the spatio-temporal power spectrum 
of wave amplitude. When wave turbulence regime is reached, we observe power-law spectra both 
in frequency and in wave number, whose exponents are found in agreement with the predictions of 
capillary wave turbulence theory. Finally, the temporal dynamics of the spatial energy spectrum 
highlights the occurrence of stochastic bursts transferring wave energy through the spatial scales. 

PACS numbers: 47.27.-i, 05.45.-a, 47.35.-i, 47.52.+j 



I. INTRODUCTION 

Wave turbulence concerns the study of the statistical 
properties of a set of numerous waves in nonlinear in- 
teraction. For strong enough interactions, a turbulent 
cascade transfers wave energy from an injection scale to- 
wards a dissipation scale. This phenomenon is analyti- 
cally described in a weakly nonlinear regime by the weak 
turbulence theory. Analytical expressions of the wave 
spectrum as a scale power-law is then derived in out- 
of-equilibrium situations in nearly all fields of physics 
involving waves (for reviews see pi, l2l). Several recent 
studies have tested the relevance of wave turbulence the- 
ory in well controlled laboratory experiments notably for 
bending waves in elastic plates [3 and for hydrodynamic 
waves on the surface of a fluid in the gravity and capillary 
regimes (for a review see [4 ). Most in situ or labora- 
tory measurements on wave turbulence involve time sig- 
nals at a fixed location and show partial agreement with 
the theory. Spatio-temporal measurements of the turbu- 
lent wave amplitude are thus needed to investigate basic 
mechanisms of wave turbulence. For gravity wave turbu- 
lence, laboratory experiments involve measurements re- 
solved in time and restricted to ID space [5 , or simulta- 
neously resolved in time and in 2D space ^ {J^ . For the 
capillary regime, most of previous laboratory works have 
tried to isolate capillary wave turbulence from the grav- 
ity wave regime by using a parametric forcing j8HT3] , by 
operating under microgravity [14 or by studying waves 
at the interface between two fluids of same density [15 . 

To our knowledge, there is no experiment studying 
capillary wave turbulence with a wave amplitude mea- 
surement simultaneously resolved in time and 2D space. 
As a consequence, linear and non linear dispersion re- 
lations has not been measured for capillary waves in a 
turbulent regime. Accurate characterization of capillary 
wave turbulence is of primmest interest at fundamental 
level. According weak turbulence theory, energy transfer 
through the scales occurs by three- wave interactions, in- 
stead four- wave interactions for gravity wave turbulence. 
Therefore, capillary and gravity wave turbulence differ 
fundamentally. Moreover, the validity domain of wave 



turbulence theory for capillary waves in experiments re- 
mains an open question, due to restrictive hypotheses of 
the theory (infinite isotropic and homogeneous system, 
weak non-linearity ...). In particular for scales smaller 
than or equal to the capillary length, viscous damping 
of waves is not negligible [16 , whereas theory lies on 
an Hamiltonian structure of the waves field. Moreover 
in oceanography, capillary waves dynamics could control 
heat and gaseous exchanges between ocean and atmo- 
sphere and contribute to the overall dissipation of gravity 
waves. 

In this article, we report full space and time resolved 
experiments on the statistics of capillary wave turbu- 
lence generated by gravity waves, a forcing close to the 
oceanographic case. The free-surface elevation is directly 
measured by a sensitive optical method known as the 
Diffusing Light Photography (DLP) [8, 9 . This tech- 
nique has a better spatial resolution and a higher sensi- 
tivity than the Fourier Transform Profilometry (an opti- 
cal technique generally used to study wave turbulence at 
the gravity wave scales) [6l|7]), and is not limited to small 
wave steepness contrary to other optical methods mea- 
suring the local wave slope [171 HH] • The DLP method 
was introduced more than 15 years ago to study capillary 
wave turbulence excited parametrically [8^ , and was able 
to get a 2D spatial measurement of the free surface am- 
plitude but only at a given time (a CCD sensor recording 
only a photography of the free surface). In our work, by 
associating DLP technique with a fast camera nowadays 
available, we present for the first time the 2D spatial 
and temporal statistics of capillary wave turbulence ex- 
cited by gravity waves. Note that recently, dynamics of 
highly non linear Faraday waves were studied using the 
DLP technique with a fast camera, but wave turbulence 
regimes have not been investigated [19] . 



II. EXPERIMENTAL SETUP 

The experimental setup is displayed in Fig. [ija). A 
Plexiglass tank (165 x 165 mm^) is filled with a diffus- 
ing liquid (IL of distilled water with SmL of Intralipid 
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20%) up to a height ho = 30 mm. Intrahpid 20 % 
(Fresenius Kabi is a commercial hpidic emulsion 
of microspheres, whose aqueous solutions are used as 
model diffusing media with characterized optical prop- 
erties [20 . Low dilution does not significantly modify 
the fluid viscosity and density from the pure water val- 
ues {u = 10~^m/s^ and p = 1000 kg/m^, respectively). 
The value of surface tension 7 is obtained from the spatio- 
temporal measurements and is found to be 7 = 60 mN/m 
(see below). The transition between gravity and capillary 
waves is expected to occur for a critical wave number 
= ypg/l which corresponds to a critical frequency 
fc ~ 14 Hz. Surface waves are generated by the hori- 
zontal motion of a rectangular paddle (130 mm in width 
and 13 mm in immersed depth) driven by an electromag- 
netic shaker (LDS V406) subjected to a random forcing 
(in phase and amplitude) band-pass filtered in frequency 
between 4 and 6 Hz. By enhancing the initial mixing of 
waves, this type of forcing is known to produce cascades 
of gravito- capillary wave turbulence in laboratory experi- 
ments [6l[Zl[22]. A LED device Phlox (lOOx lOOmm^) en- 
sures a homogeneous lighting below the transparent tank. 
A fast camera (PCO EDGE), one meter above, is focused 
on the liquid free-surface and records with 1024 x 1120 
pixels and a 200 Hz frame rate on an observation area S 
of 89 X 96 mm^. Complementary measurements are also 
performed using a faster camera (Phantom V9) with a 
1 kHz frame rate, but with a lower sensitivity, and a lim- 
ited number of images. For comparison, the wave am- 
plitude at some location is also recorded by a capacitive 
wave gauge lying outside the camera field. 

Typical 2D spatio-temporal wave elevation is measured 
as follows. If the free mean path of light within the liquid 
is smaller than the depth of the fluid and larger than the 
wave amplitude, light going through the diffusing liquid 
is no more propagating as a ray, due to the scattering 
on micrometric particles. As a result, focusing of rays 
producing caustics should disappear [9 . If the free mean 
path is also larger than the wave amplitude, the light am- 
plitude /(x, y) recorded by the camera on the free surface 
is controlled by the local depth of fluid h{x^y) [9 . After 
static calibration on a flat surface, we find in agreement 
with Xia et al. [19], an exponential decay of light ampli- 
tude with the local elevation 



/(x, y) = /o(x, y) exp {-h{x, y)/C*) , 



(1) 



Jo being the light amplitude for h = ho, and £* 25 mm 
a typical distance of order of the free mean path of light 
within the liquid. being not too small compared to the 
fluid depth ho = 30 mm, the multiple diffusion regime is 
not reached which guarantees the locality of the measure- 
ment and implies an exponential decay of light through 
the fluid [21 . The reconstruction of the free-surface to- 
pography h(x, y, t) is then performed at each time step 
using Eq. ([l]). Note that images are spatially filtered 
with a Gaussian kernel of 5 pixel size (due to the cam- 
era photon shot noise) leading to a spatial resolution of 
roughly 0.5 mm. An example of free-surface reconstruc- 
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FIG. 1. (color online) (a) Experimental setup (see text), (b) 
Snapshot of the wave field h{x,y). Random forcing between 
4 and GHz. an = 3.4mm and as = 0.30 (see text). Colorscale 
is in mm. Wave maker is parallel to the y-axis and located 
at X = —12 mm. (c) Snapshot of the spatial gradient of wave 
elevation 1 1 V/i(x, ^) 1 1 . Same parameters. 



tion is shown in Fig. [TJb). Gravity waves mainly prop- 
agate in the direction of wave maker vibrations. The 
corresponding spatial gradient of wave elevation ||V/i||, 
or wave steepness, is also shown in Fig. [ijc). Movies of 
h{x,y,t) and \\Vh\\{x,y,t) are also shown in [23]. Al- 
though the wave field is found to be inhomogeneous and 
anisotropic, capillary waves will be found in a wave tur- 
bulence regime (see below). The full space and time re- 
solved power spectrum of wave elevation Shioj, k) is com- 
puted from the set of free-surface images h{x,y,t), by 
performing successively a two-dimensional Fourier trans- 
form in space, and a Fourier transform in time, then inte- 
grating over the different directions of k. This operation 
is performed on 8192 (resp. 4196) images corresponding 
to a duration of 41s (resp. 4.2 s) for the 200 fps (resp. 
1000 fps) camera. Two important parameters are the 
typical rms wave amplitude 
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^ J h'^{x,y,t)dxdy - (^^ J h{x,y,t)dxdy ^ 



and the typical rms wave steepness 



^ JjVh{x,y,t)\\^dxdy- (^^ ljVh\\dxdy^ 



where (•) denotes a temporal averaging and a spatial 
integration on the surface S. 
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FIG. 2. (color online) Spatio-temporal spectrum of wave el- 
evation Sh{oj^k) in sinusoidal forcing with an excitation fre- 
quency of 40 Hz. Amplitude (th = 0.22 mm and Gs = 0.16. 
(+) Experimental dispersion relation extracted from the max- 
ima of 5*^(0;, /c). Red dashed line: linear theoretical dispersion 
relation. Colorscale corresponds to log^^Q (5*^ (cj, /c)). Camera 
frame rate: 200 Hz. 



III. LINEAR AND NONLINEAR DISPERSION 
RELATIONS 

Let us first focus on the linear and nonlinear disper- 
sion relation of capillary waves. For a random or a 
purely sinusoidal forcing at weak amplitude, the wave 
field should follow the linear gravity-capillary disper- 
sion relation. In order to test this assumption, spatio- 
temporal power spectrum of wave elevation Sh{(^^k) is 
displayed in Fig. |2] for a sinusoidal forcing with an exci- 
tation frequency of 40 Hz. The maxima (black crosses) 
of the spectrum Sh{oJ^k) give the experimental disper- 
sion relation corresponding to the localization of wave 
energy in the {u = 27r/, k = 27r/A) space. A is the wave- 
length and / the frequency of waves. The experimental 
dispersion relation can be accurately fitted by the linear 
dispersion relation where the surface tension 7 is the only 
free parameter 



^k'^ + (7/p) k^] tanh(/c/io) 



(2) 



g = 9.81m/s^ being the gravity acceleration and p = 
1000 kg/m^ the fluid density. We find jfu = 59.6 mN/m 
and then we consider 7 to be equal to 60 mN/m. Note 
that such a linear dispersion relation is observable, simul- 
taneously at all scales, for strong enough forcing ampli- 
tude where nonlinearity spreads energy continuously in 
the (cj, k) space. Moreover, harmonic and subharmonic 
peaks are still visible (e.g. at / = 20 Hz due to cross- 
wave instability [24 ), and the wave turbulence regime is 
thus not reached. 

At higher amplitudes of sinusoidal or random forcing, a 
departure from the linear dispersion relation is observed 
in particular in wave turbulence regimes as depicted in 
Fig. [sj^a) for a moderate amplitude and in Fig. [sjb) for 
a stronger amplitude. The curve of the maxima of the 
spectrum (black crosses) differs significantly from the lin- 
ear dispersion relation (red dashed line). These discrep- 
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FIG. 3. (color online) (a) Spatio-temporal spectrum of wave 
elevation Sh{oj^k) for two different random forcing ampli- 
tudes: (a) Gh — 2.7mm and Gs — 0.26. (b) au — 3.6 mm and 
Gs = 0.34. (+) Experimental dispersion relation extracted 
from the maxima of Sh{oj^k). Red dashed line: linear theo- 
retical dispersion relation. White dot-dashed line: theoreti- 
cal dispersion relation with a nonlinear shift (see text). Col- 
orscales correspond to log-^Q^Shioj^k)). Camera frame rates: 
(a) IkHz and (b) 200 Hz. 



ancy is found to increase with the forcing amplitude, and 
thus suggests a possible nonlinear shift of the dispersion 
relation. Indeed, dispersion relation of sinusoidal waves 
of high amplitude are theoretically known to experience 
a nonlinear shift in the gravity (Stokes waves) [25 and 
capillary regimes [26 . Although such a nonlinear shift 
is not predicted for nonlinear interacting waves with a 
continuous spectrum, an estimated nonlinear dispersion 
relation 



uj^ = \gk[l^ [akf] 



Ik' 



-1/4N 



tanh(A:/io) 



is also plotted in Fig. [3] (white dot-dashed line), a being 
the wave amplitude assumed to be equal to (Jh- The first 
term on the right-hand side of the equation (nonlinear 
gravity wave) dominates at low wave numbers and leads 
to a theoretical nonlinear dispersion relation below the 
linear one in the (cj, k) space. Experimental results fol- 
low a similar trend and suggest that nonlinear effects (ap- 
pearance of Stokes waves) are the cause of the observed 
shift of the dispersion relation. Note that bound waves 
are not observed here contrary to larger size experiments 
studying the gravity wave turbulence regime [6 . 

Directional properties of the wave field are now investi- 
gated by plotting Sh{(^^ kx^ky) as a function of A^;"^ and 
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FIG. 4. (color online) Spectrum Sh{(^,Xx~ ,Ay~ ) for fixed 
/ = 10, 20, 40 and 80 Hz {logio colorscale) showing anisotropy 
of wave field at high forcing amplitude (an = 3.6 mm and 
as = 0.34). Dashed black circle: Linear dispersion relation of 
Eq. Q . Wave amplitudes are larger in the direction of forcing 
(x-axis) revealing anisotropy in all frequency range. 



Xy~^ for fixed / = co'/(27r) as in Fig. [ij The direction 
of the wave forcing being along the x-axis, we observe a 
strong anisotropy of the wave field. This anisotropy is 
conserved whatever the frequency scales inside the capil- 
lary waves range (see at / = 20, 40 or 80 Hz). Capillary 
waves appears preferentially along the same direction as 
the long waves from which they take origin by nonlinear 
interactions. Due to the geometry of our experimental 
device and to the viscous dissipation, wave field appears 
also inhomogeneous. The spatial inhomogeneity of the 
wave field can be estimated by taking the square root of 
the temporal average of h^{x^y^t) and then computing 
the spatial standard deviation. We found in our exper- 
iments a spatial inhomogeneity of order 24%, without 
significant changes with the wave amplitude. 



IV. TEMPORAL AND SPATIAL SPECTRA OF 
WAVE ELEVATION 

We next investigate the scaling of the spectrum with 
the spatial and temporal scales. Theoretically, the spec- 
trum of capillary wave turbulence is predicted to scale 
as S'f ^^(cj) - cj-i^/^ and as 5f ^^(A:) - A:"!^/^ [27^,^281. 
Experimentally, the temporal spectrum of wave eleva- 
tion Sh{oj) is obtained by averaging Sh{oJ^k) over k as 
shown in Fig.jsja) for different forcing amplitudes. When 
the forcing amplitude is increased, the peaks related to 
the forcing frequencies and their harmonics progressively 
disappear on the spectra, and a power-law spectrum is 
observed in the capillary regime on one decade in fre- 
quency (20 < / < 200 Hz). No self-similar regime is 
observed in the gravity range due to the small size of 
the tank. Exponents from frequency power-law fits of 
spectra in the capillary range are plotted in Fig. jsjb) as 
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FIG. 5. (color online) (a) Temporal power spectra Sh{(^) for 
different forcing amplitudes. From bottom to top: an =1.3, 
2.7, 2.1, 3.4 and 3.6 mm, and as 0.15, 0.26, 0.19, 0.3, and 
0.34. Solid black line is the capillary prediction f~^^/^ ^ The 
second measurement is recorded with a 1 kHz frame rate, the 
other ones at 200 Hz. (b) Inset: 5*^ (cj) -exponents vs. ah 
(fits from 20 < / < 100 Hz). Dashed line shows the theo- 
retical value —17/6. (c) Spatial power spectra Sh{k) for the 
same measurements. Solid black line is the capillary predic- 
tion k~'^^^^. (d) Inset: 5'^(A:)-exponents vs. an (fits from 
0.094 < < 0.30 mm~^). Dashed line shows the theoret- 
ical value —15/4. is a characteristic scale of dissipation. 
Ac is the gravity- capillary crossover and = 0.094 mm 
(See Fig[6|. (e) Spatial spectrum of the energy E{k) (•). 
Contributions of gravitational Egrav ^ Sh{k) and capillary 
Ecap ~ k^ Sh{k) energies are well separated. The solid 
black line is the prediction for the capillary energy spectrum 
^ A;~^/^. ah — 3.6mm and as — 0.34. 



a function of the rms wave amplitude, These spec- 
tral exponents are in a good agreement with the pre- 
dicted exponent f-^^l^ at high enough wave amplitude. 
Note that the signal (recorded with the highest frame 
rate) reaches noise level around 200 Hz before reaching 
the dissipative range. These results are confirmed by si- 
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multaneous measurements provided by a capacitive wave 
gauge, which was used extensively in previous studies of 
capillary wave turbulence [4| [22] . 

The spatial spectrum Sh{k) is obtained by averaging 
Sh{oJ^ k) over the frequencies /, and is shown in Fig.jsjc) 
for different forcing amplitudes. For high enough wave 
amplitude, a power-law spectrum is also observed in the 
capillary range. Exponents from wave number power- 
law fits of spectra in the capillary range are plotted in 
Fig. [5|d) and are close to the weak turbulence predic- 
tion in k~^^^^. Finally, note that the departure from 
the theoretical power-law at small scales is observed in 
Fig. [sjc) at \d~^ • y^d~^ is found to increase with an 
increasing forcing amplitude. This suggests that A^^"^ 
arises from the balance between viscous dissipation and 
nonlinear interactions [l6l|27]. Moreover, for the highest 
forcing amplitude, A^^" = 0.48 mm~^ which corresponds 
to fd ~ 200 Hz. For A~^ > A^^"^, self similarity is bro- 
ken and spectrum departs from the power-law. Although 
some strong hypotheses of weak turbulence (notably ho- 
mogeneous, isotropic and infinite system) are not verified 
in our experimental system, a correct agreement is found 
for both the temporal and spatial scalings of the capillary 
wave spectrum. 




V. STATISTICS OF ENERGY MODES 

We now use the energetic point of view to better probe 
the statistics of capillary wave turbulence. Assuming 
that, in average on a wave period, there is equality be- 
tween the kinetic and potential energies at each scale, 
the total wave energy at one scale should thus be pro- 
portional to the potential energy at this scale. There- 
fore, neglecting non linear surface deformations, one can 
compute the total potential energy of the wave field as 



-'pot 



n 



pg\h{k)\'^je\h{k)\' 



dk^ J E{k)dk. 



From the spatial power spectrum Sh{k) ~ \h{k)\'^ {h 
denotes the Fourier transform of /i), one can thus eas- 
ily deduce the energy spectrum E{k) that is depicted in 
Fig. [sje). The separation of the respective contribution 
of the gravity energy and of the capillary energy appears 
clearly. Moreover, in the capillary range, E{k) follows 
the predicted spectrum of weak turbulence for the en- 
ergy ~ confirming the above results. 

The spatial power spectrum Sh{k) is now computed at 
each time step from the wave field h{x,y^t). We thus 
obtain the temporal evolution of wave energy spectrum 
E{k^ t) as shown in Fig.[6]^a). An important result is that 
E{k^ t) exhibits stochastic bursts transferring energy to- 
wards small spatial scales. Indeed, most of the time, 
the wave energy is confined near the forcing scales, but 
from time to time, energy cascades through all spatial 
scales (within the inertial capillary range). Moreover, 
at a given spatial scale /c^, i.e. an horizontal slice in 
Fig. [6]^a), E(k^^t) is found to be strongly erratic and 
bursts of random large- amplitude occur as for gravity 
wave turbulence [6]. The probability density function 
(PDF) of E{k^^t) of a mode k^ is shown in Fig. |6|b) 
for two forcing amplitudes. For a moderate forcing, the 
PDF roughly follows an exponential distribution as ex- 
pected for a Gaussian wave field [2 . At higher forcing, 
a departure from the exponential distribution becomes 
significant, likely revealing a stochastic nonlinear phe- 
nomenon as described theoretically [29] . The PDF shape 
does not depend significantly on the scale k^ within the 
inertial range of the capillary cascade. 



VI. CONCLUSION 



FIG. 6. (color online) (a) Evolution of the wave energy E{k^ t) 
in the (t,A~^) space. Color scale in log^^Q scale, (b) Probability 
density function of wave energy E(k^, t)/ {E(ki,, t)) at a mode 
— 0.094 mm~"^: (*) moderate forcing (an — 2.0 mm, as — 
0.19), and (o) high forcing (an = 3.6 mm, as = 0.34). Dashed 
line is the exponential distribution: exp(— x). 



We have reported experiments on statistics of capillary 
wave turbulence generated by gravity waves. The full 
space and time resolved spectrum of wave height is 
obtained for the first time in the capillary range by 
means of an optical method. The spectrum exhibits both 
frequency and wave number power-law scalings whose 
exponents are close to the ones of the weak turbulence 
theory. Beyond the confirmation of these predictions 
in the capillary regime, significant departure from the- 
oretical hypotheses has been reported (notably related 
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to the isotropy and homogeneity of the wave field). We 
have also observed the occurrence of stochastic bursts in 
time transferring wave energy through the spatial scales 
within all the inert ial range. The joint space and time 
resolved measurement have also allowed us to observe 
the linear and non linear dispersion relations of waves in 
a turbulent regime. When the forcing is strong enough, 
a nonlinear shift of the dispersion relation is reported, 
the spectrum power-law scalings with k and uj being still 
in agreement with the predictions. Although the latter 
are derived using the linear dispersion, nonlinearity 
is considered in the theory as a perturbation of the 
linear state, and the observations reported here are thus 
consistent with the predictions. Further studies will 
investigate more in detail the statistics of Fourier modes 



of wave energy as well as the space-time correlations in 
order to better understand basic mechanisms involved 
in capillary wave turbulence. 
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